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Abstract 

A mathematical method for through-wall imaging via wave phenomena in the time do¬ 
main is introduced. The method makes use of a single reflected wave over a finite time 
interval and gives us a criterion whether a penetrable obstacle exists or not in a general 
rough background medium. Moreover, if the obstacle exists, the lower and upper estimates 
of the distance between the obstacle and the center point of the support of the initial data 
are given. As an evidence of the potential of the method two applications are also given. 
AMS: 35R30 
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1 Introduction 

Assume that there is a large wall between an observer and an unknown object. The wall is opaque 
and thus the observer can not see the object directly. How can the observer find the object? 
Consider the case when the wall is electromagentically penetrable. Then there is no dought that 
everyone thinks about the use of electromagnetic wave. Generate the electromagnetic wave at 
the place where the observer is and observe the reflected wave from the object at the same place 
over a finite time interval. The observed wave should include information about the object. 
How can one extract the infromation from the observed wave? The prurpose of this paper is 
to develop a mathematical method for throuh-wall imaging by using the governing equation of 
the wave from the begining to end. Originally the governing equation should be the Maxwell 
system. In this paper, as a first step we consider two single equations for scalar waves. 


1.1 Finding discontinnity in a medium with a rough refractive index 


Let us formulate the first problem. In this paper we consider only the case when re = 1,3. Let 
a € L°°(R”) and satisfies ess.inf2,gR,nQ;(x) > 0. Given / G L^(R"’), let u = Uf{x,t) be the weak 
solution of of 


a{x)d^u — Are = 0 inR^'x ]0, T[, 


( 1 . 1 ) 


re(x,0) = 0 inR"’, dtu{x,0) = f{x) inR"’. 


The notion of the weak solution is taken from [3]. See Subsection 2.1 in this paper. 
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We assume that a takes the form 


a{x) 


ao{x), if X G R"" \ D 
ao{x) + h{x), X G D, 


( 1 . 2 ) 


where D C R” is a bounded open subset with Lipschitz boundary. The function ao belongs to 
L°°(R"') and satisfies ttiq < ao{x) < Mq a.e .X G R” for positive constants mo and Mq] h on D 
belongs to L°°{D) and satisfies one of (A.I) and (A.II) listed below: 

(A.I) BC > 0 h{x) > C a.e. x G D] 

(A.II) 3C > 0 — h{x) > C a.e. x G D. 

D is a model of an unknown penetrable obstacle with material parameter aQ + h embedded 
in the background medium with material parameter aQ. The distribution of ccq models various 
penetrable walls in the space and D is something hidden in the walls or a space between the 
walls and various penetrable obstacles. 

Let B be an open ball satifying BnD = ij). Let / G L^(R”) satisfy f{x) = 0 a.e. x G R^ \ B 
and ess.infajgs/(x) > 0. Generate w/ by the initial data /. In this paper, we consider the 
following inverse problems under the assumption that oq is known and that both D and h are 
unknown. 

Problem 1. Find a criterion whether = 0 or not in terms of only Uf on B over time interval 
10. T[. 


Problem 2. Assume that D Extract information about D from Uf on B over time interval 
]0, T[ for a fixed large T. 

Note that the correspondence {D,h) i— u\bx]o,t[ is nonlinear and thus both problems 
become nonlinear problems. The existence of variation of oq outside D plays a role of the 
wall in which a wave can propagate. This paper aims at developing a through-the wall imaging 
method via the governing equation of the wave used. See, e.g., [1, 19] and references therein for 
practical approaches. 

In this paper, we employ the Enclosure Method introduced by the author himself in [6, 7] 
as a guiding principle for attacking the problem mentinoned above. It is a direct method which 
connectes the unknown discontinuity and the observaton data through the so-called indicator 
function. In [8] the author initiated to apply the method to inverse obstacle problems whose 
governing equations are given by the heat and wave equations in one-space dimension. Now we 
have several applications of the method to inverse obstacle scattering problems in three-space 
dimensions whose governing equations are given by the wave equautions [9, 10, 11, 12]. See [13] 
for a review of the results in those papers. However, imaging an obstacle through-the wall is 
still out of the range of the results in those papers. 

Now let us describe the result. 

Let r > 0 and define 



e '^^u{x,t)dt inR"". 


(1.3) 


Let V = v{x,t) G iL^(R"') be the weak solution of 


Av — aQv'^v + aof = 0 inR”. 


(1.4) 


The following resut gives us solutions to Problems 1 and 2. 

Theorem 1.1. We have: 

(i) if D = (h, then for all T > 0 it holds that 

lim / aQfiw — v)dx = 0; 

r^oo Jb 
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(ii) if D ^ (b and h satisfies (A.I), then for all T > 2Modist{D, B) it holds that 

lim e'^'^ / aof(w — v)dx = —oo: 

r-S-oo 

(Hi) if D 9 and h satisfies (A.II), then for all T > 2MQdist{D, B) it holds that 


lim e'^'^ / aof(w — v)dx = oo. 

— Jb 


Moreover, in case of both (ii) and (Hi) we have, for all T > 2MQdist{D, B) 


and 


lim inf — log 

T ->00 2t 


lim sup — log 

T->00 2t 


IB 


aof{w — v)dx 


IB 


aof{w — v)dx 


> —MQdist{D, B) 

< —modist {D, B). 


(1.5) 

( 1 . 6 ) 


Note that ao and h are just essentially bounded on and D, respectively. In particular, 
ao can be a model for various background media such as multilayered media with complicated 
interfaces or unions of various domains with different refractive indexes. It will be impossible 
to apply any approach based on the idea of geometrical optics to this case. See [5, 15, 16, 17] 
for such approach in the case when the scattering kernel which is the observation data in the 
Lax-Phillips scattering theory is given under the assumptions that Oio{x) = 1 a.e. x G dD 
is smooth and h G C°°{D). 

Let p be the centre of B and rj the radius. Since dist {D, B) = dgnip) — estimates (1.5) and 
(1.6) give us an upper and lower estimate of dgnip) provided mo and Mq are known and T is 
sufficiently large. Thus we obtain the upper and lower estimation of the sphere \x — p\ = dQuiv) 
whose exterior encloses the object. Estimates (1.5) and (1.6) suggest a new direction of the 
Enclosure Method in the case when the background medium is inhomomegeneous and quite 
complicated; give up to find a precise quantity in the observation data which is related to the 
exact location of unknown obstacles; instead give lower and upper estimates rigorously like (1.5) 
and (1.6) for dist {D,B). 

Some further remarks are in order. 

• In Theorem 1.1 it suffices to know u on R not whole v. However, without knowing ao 
outside B it is impoosible to compute u on R in advance. In the last section we suggest an 
experimental computation procedure of u on R without detailed knowledge of ao outside R. 
It seems that this will be useful for the daily security of the interior of the room, house and 
building which have complicated environment. For the purpose Theorem 1.1 will be suitable 
since ao is just essentially bounded and we never assume any other regularity. 

• If ao(a;) = 1 a.e. x G R”, then one can choose Mq = mo = 1 and thus (1.5) and (1.6) 
imply that 

f{w — v)dx 

This coincides with a result (1.20) in [10]. 

In the next subsection we apply the idea developed here to more special case and show (1.5) 
and (1.6) can be replaced with a single formula. 


= —dist (R, R). 


lim — log 
-—>oo 2r 


IB 
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1.2 Finding discontinnity in a dissipative medium 

Let q G satisfy q{x) > 0 a.e.x G R"". Given / G L^(R"'), let u = u{x,t) be the weak 

solution of 

d^u — Au + q{x)dtu = 0 inR^x ]0, T[, 

(1.7) 

u{x,0) = 0 inR”, dtu{x,0) = f{x) inR"'. 

We assume that q takes the form 

f if X G R^ \ H 

q{x) = i (1.8) 

[ qo{x) + h{x), if X G H, 

where H is a bounded open set of R” with Lipschitz boundary. The function go belongs to 
L°°(R"') and satisfies qo{x) > 0 a.e. x G R"". 

The function h on D in (1.8) has to belong to G L°°{D). We assume that h satisfies one of 
(A.I) and (A.II). 

Let r > 0. Define w by the formula (1.3) in which u is replaced with the solution of (1.7). 
Let V G 77^ (R”) be the weak solution of 

(A — — Tqo)v + / = 0 inR”. (1.9) 


The following result is new and suggests that, even in a rough dissipative medium the solution 
of (1.7) carries information about D clearly like the wave equation outside D. 

Theorem 1.2. Let B and f be the same as those of Theorem 1.1. We have: 

(i) if D = ^, then for all T > 0 it holds that 

lim e'^'^ / f{w — v)dx = 0; 

T ^CxD J Q 

(ii) if D ^ ^ and h satisfies (A.I), then for all T > 2dist{D, B) it holds that 

lim e^'^ / f(w — v)dx = —oo: 

T—>oo 

(Hi) if D ^ (b and h satisfies (A.II), then for all T > 2dist{D,B) it holds that 

lim / f(w — v)dx = oo. 

r^oo 

Moreover, in ease of both (ii) and (Hi) we have, for all T > 2dist{D,B) 


lim — log 

T->00 2t 



v)dx 


— dist {D, B). 


( 1 . 10 ) 


Roughly speaking, we see that (1.9) as r —)• oo corresponds to (1.4) with uq = 1 and thus 
formula (1.10) is reasonable. 

Comparing Theorem 1.2 with the previous results in [10, 11] for the wave equation outside 
D, we see (ii) and (iii) suggest us that that assumptions (A.I) (stronger dissipation) and (A.II) 
(weaker dissipation) correspond to the Dirichlet and Neumann boundary conditions imposed on 
dD, respectively. 

A brief outline of this paper is as follows. In Section 2 first we specify the meaning of the 
weak solution used in the formulation of the problems. Second we establish two basic integral 
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identities for w given by (1.3). Theorem 1.1 is proved in Section 3. The proof starts with 
deriving the lower and upper estimates for the indicator function 


T 



v)dx 


as r —7- oo from the basic identities in Section 2. Next, by virtue of the governing equation 
(1.4) we see that v has pointwise explicit lower and upper esimates. Combining those, we obtain 
all the conclusions stated in Theorem 1.1. The proof of Theorem 1.2 is given in Section 4 which 
is a combination of a reduction and similar argument done in the proof of Theorem 1.1. In 
Section 5 we give some remarks concerned with a “practical use” of Theorem 1.1 and present a 
sharp form of Theorem 1.1 in the case when the space dimension is 1 and ao has a special but 
important form. 


2 Preliminaries 


2.1 On the solution class 


In this subsection we specify the meaning of the weak solutions of (1.1) and (1.7) at the same 
time. It is based on Theorem 1 given on p.558 in [3] which we have already used for the 
formulation of the weak solution of (1.1) in [10]. 

Set V = and H = L^(R"'). Applying Theorem 1, we know that given vP and 

v} G H, there exists a unique u satisfying 

u€L\0, T-,V), u' = ^€L\0, T-,V), A(C'(u'( •)) G i'(0, T;V'), 
such that, for all (/> G 14 

< ■^C{u'{t)),(J) >+bo{u'{t),v) + a{u{t),(p) = 0, a.e. tGjO, T[, (2.1) 

and u(0) = and u'(0) = u^, where 


a{u, v) 



Vtt • Vvdx, u,v £ V, 


bo{u,v) = / quvdx, u,v£V, 

JR" 

and C : H —)■ H is the bounded linear operator defined by 

C{u) = au, u £ H. 

Note that this C satisfies (5.11) on p. 553 in [3] under the condition ess.inf 2 ;gRna(x) > 0. Since 
q > 0, bo satisfies (5.8) on p. 553 with /3o = 0 in their notation. However, equation (2.1) is 
homogeneous, i.e., the source term is 0, and by virtue of this, their proof also covers this case. 

In this paper, we say that this u for = 0 and = f with q = 0 and a = 1 is the weak 
solutions of (1.1) and (1.7), respectively. 

We see that w given by 


w= e 
Jo 


^—Tt 


udt 


belongs to V and applying integration by parts to (2.1) multiplied by e we obtain, for all 
<j)£V 

I {Vw ■ V4> + {r'^a + Tq)w(j)}dx — [ af<j)dx + e~'^'^ [ J'(f>dx = 0, 
jR" JR" JR3 
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where 


-F(x) = a{x)u'{x, T) + {a{x)T + q{x))u{x, T). 
This means that w is the weak solution of the following equation: 

(A — — qT)w + af = F inR*^. 

2.2 Two basic identities 


( 2 . 2 ) 


/ /{(«o ~ oi)v + aR}dx = T^ (oq ~ a)v^dx 

J'RP JR” 

+ [ (|Vi?|^ + + 6“”^^ ( [ aFRdx— [ aFvdx] , 

jR” \JR" jR” / 

where R = w — v. 

Proof. From (1.4) We have 

/ aofwdx = T^ aQVwdx + / Vv ■ Vwdx. 

Jr" Jr” Jr” 

On the other hand, from (2.3) we have 

/ afvdx = T‘^ / awvdx + / Vw ■ Vvdx + e~'^'^ / aFvdx. 

JR” JR” JR” JR” 

Therefore we obtain 

/ f{aQW — av)dx = / (oq ~ a)vwdx — e~'^'^ / aFvdx. 

jR” jR” JR” 


Write 


/ (ckq ~ a)vwdx = / (oo ~ a)v‘^dx + / (oo ~ a)vRdx 

JR” JR” JR” 


/R” jR” 

It follows from (2.3) and (1.4) that R satisfies 

Ai? — ar^R + (oq — a)T‘^v — (oq — «)/ = ae~'^'^F inR" 


(2.3) 


In this subsection we consider only the case when q{x) = 0 a.e. x G R”. Then, it follows from 
(2.2) that w satisfies 

Arc — + af = ae~'^'^F inR”, 

where F{x, r) = v!{x, T) + rtt(x, T). 

In what follows / just belongs to L^(R"'); v G R^(R”^) satisfies (1.4). 

Proposition 2.1. We have 


(2.4) 


(2.5) 

( 2 . 6 ) 

(2.7) 


Thus we have 

[ (ao — a)vRdx = [ (|Vii|^ + + [ {ao — a)fRdx + [ aFRdx. 

jR” jR” jR” jR” 

Subsituting this into (2.6) we have 

/ {a^ — a)vwdx = / (oq — a)u^(ix + / (|Vi?p + ar^i?^)(ix 

JR” JR” JR” 

+ / (oo ~ CK)/-Rc?a^ + / aFRdx. 

JR” JR” 
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A combination of this and (2.5) yields 

/ f{aow — av)dx = T'^ {aQ — a)v'^ dx + / {\V R\'^ + R'^)dx 

Jr" Jr" Jr" 

+ / («o ~ a)fRdx + ( / aFRdx — / aFvdx ) . 

JR" VJR" JR" / 

Since {aow — av) — (ao ~ o)-R = (oo ~ oi)v + aR, we obtain (2.4). 

□ 

Proposition 2.2. VFe have 

/ /{(a — ao)'*^’~ / — (a — 

jR" jR" a 

+ J ^|Vi?p + ar^ i? + ^ n dx + e~'^^ aFRdx + j aFvdx"^ , 

where R = w — v. 

Proof. From (2.3) we have 

/ afvdx = / awvdx + / Vrc • Vvdx + / aFvdx. 

JR" JR" JR" JR" 

On the other hand, from (1.4) we have 

/ aofwdx = / a^vwdx + / Vn • Vwdx. 

JR" jR" jR" 

Therefore we obtain 

/ f{av — aow)dx = / (a — ao)nt(;(ix + / aFvdx. 

jR" jR" jR" 

Write 

/ (a — ao)vwdx = / (a — / (ao ~ a)wRdx. 

JR" JR" JR" 

From (2.7) we see that i? satisfies 

AT? — + (ao — a)T‘^w — (oq — a)/ = ae~'^'^F inR"". 

Thus we have 

/ {ao — a)wRdx = [ {\VRl"^ + aor'^R^)dx + [ {ao — a)fRdx + e~'^'^ [ aFRdx 

jR" jR" jR" jR" 

Subsituting this into (2.10) we have 

/ {a — ao)vwdx = / (a — + / (|VRp + 

jR" jR" jR" 

+ / {ao — a)fRdx + e~^'^ / aFRdx. 

JR" JR" 

Now a combination of this and (2.9) yields 

/ f{av — aow)dx = j (a — ao)w‘^dx + / (Vi?|^ + aoT‘^R?)dx 

jR" jR" jR" 


( 2 . 8 ) 


(2.9) 

( 2 . 10 ) 


+ / {ao — a)fRdx + e / aFRdx + / aFvdxj. 
JR" VJR" JR" / 


( 2 . 11 ) 
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Finally write {av — a^w) — (oq ~ oi)R = (a — ao)w — a^R and 


+ {a — aQ)w‘^ = a 

Combining these with (2.11) we obtain (2.8). 

□ 


ao 

w - V 

a 


H- [a — ao)v . 

a 


3 Proof of Theorem 1.1. 

First we derive two asymptotic estimates. 
Proposition 3.1. We have, as r —)• oo 


[ faQ{w — v)dx<T^ [ —{ao — a)v‘^dx + 0{T ^e (3.1) 


and 


/ fao{w — v)dx > / {aQ — a)v‘^dx + 0 {t ^e 

Jrp JR" 


(3.2) 


Proof. We start with describing two simple estimates for v and R = w — v. It follows from (1.4) 
that 

I {\Vv\‘^ + — aQfv)dx = 0, 

jR" 

that is, 

L {l Vop + a„ (™ - ^) I = jL 

It is easy to see that from this we obtain, as r —)■ oo 

[ {\Vv\‘^ + aoT‘^v‘^)dx = 

jR" 

and in particular. 


^IIl2(r») - 0{t ^). 


(3.3) 


Next rewrite (2.4) as 


[ {ao — a)v‘^dx+ [ \\V R]"^ + a i tR— 

jR" jR" ( \ 


/ - e--^F 
2 t 


dx 


= / (oo — a)fvdx + —j / a{f — e F)‘^dx. 
jR" 4r^ jR" 

Since F = u'{x, T) + tu{x, T) and ||u'( •, T)\\i 2 (^n-^ + ||rt( •, F)|h 2 (Rn) < oo, it follows from this 
and (3.3) that 


J ||Vi?p + a ^rii- 


/ - 

2^ 


dx = 0 {t "^) 


and hence 


[ {\VRf + aT‘^R^)dx = 0 {t ^). 

jR" 


In particular, we have 


l|-R||L2(Rn) — 0(r 


-2^ 


(3.4) 



Applying (3.3) and (3.4) to (2.4) and (2.8), we obtain, as r —)■ oo 


Here note that 


/ /{(«o — Q;)^^ + aR}dx = T^ (ao — a)v^dx 

JR" JR" 

+ / {\VR\^+ aT‘^R^)dx + 0{T-^e-^^) 

JR" 

/ /{(a — ao)^'^ “ / — (a — ao)'f^^'ia^ 

JR" JR" a 

+ J ^|Vi?p + ar^ i? + ^ V ^ dx + 0(r“^e“'^^). 


, . ao. \ ~ “o)^ 

(^CKq — Qij “h -— ^0/ —- 

a a 


(ao — a)r; + aR + (a — aQ)w — a^R = 2(a — ao)i?. 
Thus, summing (3.5) and (3.6) up, we obtain 

2 [ f(a — an)Rdx + f - - 


2 [ /(a — ao)Rdx + f -— v'^dx 

JR" JR" a 

= [ (|Vd?p + ar^i?^)dx + [ [|Vi?p + ar^ ii+(^l ——) dx + 0(r“^e“'^^). 
JR" JR" \ \ a J J 


Now assume that a = ao on supp /. Then (3.8) gives 


HVRl"^+ aT‘^R^)dx + |^|Vi?|2+ ar2 i?+(^1 - X ^ dx 

2 f (a — ao)^ I -1 -tT\ 

= T / -X dx + U{t e ) 

jR" a 


and hence 


Note that 


[ (|Vi?p + aT^ii^)dx < [ 

jR" JE 


2/ («-«0)% 2^x + O(T-^e-^). 

JR" a 


(a - ao) 


-h (ao — a) = —(ao — a). 

a a 

Then a combination of (3.9) and (3.5) gives (3.1) and (3.2), respectively. 

□ 

Remark 3.1. We note that (3.7) gives 

ao — a < — (ao — a), 
a 

Thus (3.1) and (3.2) are reasonable. 

Now assume that iJ = 0. Then ao = a and thus (3.1) and (3.2) yield 


[ faoRdx = 0 {t '^^). 

JR" 
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This gives (i). 

The proof of (ii) is as follows. Since h satisfies (A.I), it follows from (3.1) and (3.2) that 


j^^faoRdx <-At‘^\\v\\\2(^i:))+0{t ^'^), 

where A = —(mg)/(M^)C', M = ^ ess.sup 2 ,gj^na(x) and C comes from (A.I); 

faoRdx > -A'r^||x||| 2 (j)) + 0(r"^e"^^), 


(3.10) 


(3.11) 


where A! = 


L°°(D) 


> 0 . 


Thus it suffices to give a lower and upper estimate of v over D. Given A > 0 define 


1 


KxiO = 



if n = 1, 

e-AkI 
47r|^| ’ 

if n = 3. 


It seems that the following lemma is closely related to the maxiumu principle or comparison 
principle [4]. However, our final purpose is to consider the electromagnetic wave which satisfies 
a system. So in Appendix we give a proof without making use of such principles. 

Lemma 3.1. Let f G L^(R'^) and satisfy f{x) > 0 a.e. x G R". Let v G R^(R") be the weak 
solution of (1-4)■ We have 


and 


v(x) > / ao{y)f{y)KMQT{x - y)dy a.e.x G R'", 
JR" 

v{x) < / ao{y)f{y)KmoT{x - y)dy a.e.x £ R*". 
JR" 


(3.12) 

(3.13) 


Let us continue the proof of (ii). Since the case when n = 1 can be treated 
we only consider the case when n = 3. By the mean value theorem [2] we have 


easily, hereafter 


1 r (firMor]) e-^^o\x-p\ 

An Jb \x — y\ ^ {Mqt)^ \x — p\ 


(3.14) 


where p and rj are the centre and radius of B, respectively and = ^cosh^ — sinh,^. 

By [14], we know that there exists a positive constant C and number G R such that, for 
all r >> 1 

r p~2t|x— p| 

1 > c'. 

Jd \x -p\^ 

Applying this and (3.14) to a lower bound derived from the lower bounds for ao and f on B for 
the right-hand side on (3.12), we obtain, for all r >> 1 


(3.15) 

2wmV(iWo^ 2MorJau(p) 

- (Mor)6+M 
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where C = ess-inf^^gs f{x)- Since as ^ 


oo 


from (3.15) we obtain, for all r >> 1 

where C" is a positive constant. Since dQj:,{p) — r] = dist {D, B) we finally obtain, for all r >> 1 


> 


\M\l^(d) — 

Now a combination of this and (3.10) yields 


(^W.j.-(4+At)g-2Mordist (D,B) 


(3.16) 


lim / faoRdx = —oo 

-- >00 J^n 


if T > 2Modist {D, B). 

Moreover, it is easy to obtain from (3.13) that, for a positive constant C 
have 

^ -2moTdist (D,B) 


and all r > 0 we 
(3.17) 


Now it is easy to verify that a combination of (3.10) and (3.16) yields (1.5); a combination of 
(3.11) and (3.17) yields (1.6). 

A similar argument based on (3.1) and (3.2) works also for the case when h satisfies (A.II). 


4 Proof of Theorem 1.2 

From (2.3) we have 


where 

Define 


(A - - Tq)w + f = e-^^F{x, r) in R" 

F{x, t) = u'{x, T) + (r + q{x))u{x, T). 
a = a{x, r) = 1 + —, 


f = fix,T)= 1 + ^ / 


and 


Then (4.1) becomes 


-1 


-1 


F = F(x,t) = ( 1 + T ) F. 


(A — T‘^a)w + af = de F inR”. 


Let V G R^(R"') be the solution of (1.9). Define 


cfo = C(o{x,t) = 1 + — 
r 


and 


/o = fo{x,T) = IH-/. 

r 


qo 


-1 


(4.1) 
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Then (1.9) becomes 


(4.2) 


(A — doT‘^)v + dofo = 0 mR”. 

Note that ao/o = / = «/ and dF = F. 

Define R = w — v. We start with the following two integral identities which can be derived 
along the same line as that of Propositions 2.2 and 2.3 and thus we omit the description of their 
proofs: 

[ fRdx = [ {do — d)v^dx + [ (|Vi?p + dT‘^R^)dx 

JR" JR" JR" 


+e 


-tT 


FRdx — 


Fvdx 


'R" 


/R" 


/R’ 


fRdx = T^ [ —{d — do)v‘^dx 

JR" a 


+ MVRp + dr^ 


/R" 

Note also that 


and 


R+ ( 1-^ It; 

a 


dx + e 


-tT 


FRdx + / Fvdx 


'R" 


/R" 


/ (cfo ~ d)v‘^dx = T {qo — q)v‘^dx 

JR" JR" 

_^2 f ^ — j- f ''' F qo _ q)v‘^dx. 

jR" a jR" T + q 


(4.3) 


(4.4) 


/R" T + q 

Then, applying a similar argument as done in the proof of Proposition 3.1 to (4.3) and (4.4) we 
obtain 


Proposition 4.1. We have, as r —)• oo 


and 



fRdx <t[ 
JR" 



fRdx > T 


(go - q)v‘^dx + 0{t ^^) 

T + q 

f {qo - q)v'^dx + 0{T~^e~'^'^). 

R" 


(4.5) 

(4.6) 


Thus it suffices to prepare the following lower and upper estimates for v. 

Lemma 4.1. Let f G L^(R”) and satisfy f{x) > 0 a.e. x G R'^. Let v G H^{IV^) be the weak 
solution of (1.9). We have 


and 


where 


v{x) > / f{y)KL (^^).,{x - y)dy a.e. x G R'" 
JR" 

v{x) < / f{y)Kr{x - y)dy a.e. x G R'", 
JR" 


Lo{t) 



(4.7) 

(4.8) 


(4.9) 


and Lo = ess.sup^^-^n qo{x). 

Proof. Since v satisfies (4.2) and Lo{t) given by (4.9) satisfies cfo(x) < To(t)^ a.e. x G R”, from 
(3.12) in Lemma 3.1 one obtains (4.7). The proof of (4.8) is as follows, v has the expression 


V = vo- eo, 
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where vo G H^CR^) solves 

(A — t‘^)vo + / = 0 inR*^ 

and eo G i7^(R"') solves 

(A - - Tqo)eo + rq^vo = 0 InR"-. (4.10) 


vq has the explicit form 

vo{x) = [ f{y)Kr{x-y)dy>0 
JR" 

and hence rqovo > 0. Applying (4.7) to (4.10) we obtain 

eo(a:) > / Tqo{y)vo{y)KL (^^)^{x-y)dy 

JR" 

and hence €o(x) > 0 a.e. x G R”. Therefore we obtain v < vq and thus (4.8). 

□ 

Now noting Lo(t) —)• 1 as r —)■ oo, it is not difficult to deduce all the conclusions in 
Theorem 1.2 from (4.5), (4.6) and Lemma 4.1 as done in the proof of Theorem 1.1. 


5 Remarks and further problems 

5.1 How to compute the indicator function without knowledge abont oq ont- 
side B 


In this work we think that D is an invader into a space with the refractive index aQ. However, it 
will be difficult to know the detail of oq and find u on H which is the solution of equation (1.4). 
In this section we describe an experimental computation procedure of the indicator function 
from the observed data in the two spaces one of which has an invader and another does not have 
an invader yet. 

Let V = V{x,t) be the weak solution of 


aoix)d^V -AV = 0 inR^x ]0, T), 


l/(x,0) = 0 inR"", dtV{x,0) = xb{x) inR"". 


(5.1) 


Generate wave 14 by the initial data, say / = xb in the space which has no invader and observe 
Ve on B over time interval ]0, T[. We assume that 14 on H is given by V on B. 

From Proposition 3.1 in the case when ao = a we obtain 


where 


Thus, we obtain 


/ ao{ve — v)dx = 0(r 
JB 

pT 

Ve(x,T) = / e~'^'^Ve{x,t)dt. 

Jo 


/ ao{w — Ve)dx = / ao{w — v)dx + 0 {t ’^^). 

Jb Jb 

Then, it follows from Proposition 3.1 that, as r —)• oo 



Ve)dx < / — (ao 

JR" OL 


a)v^dx + 0(r ^^) 
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and 


/ ao{w — Ve)dx > / {ao — a)v‘^dx + 0 {t 

JB JR" 


Therefore, one can transplant all the results in Theorem 1.1 into the present case and we 
obtain 

Theorem 5.1. We have: 

(i) if D = ^, then for all T > 0 it holds that 

lim / aof(w — Ve)dx = 0; 

T^OO 

(a) if D ^ (b and h satisfies (A.I), then for all T > 2Modist{D, B) it holds that 

lim / aofiw — Ve)dx = —oo; 

T^OO 

(Hi) if D ^ (b and h satisfies (A.II), then for all T > 2MQdist{D, B) it holds that 

lim / aof{w — Ve)dx = oo. 

T ^OO J Q 

Moreover, in case of both (ii) and (in) we have, for all T > 2MQdist{D, B) 

o:of{w — Ve)dx >—MQdist{D, B)] 

Oiofiw — Ve)dx <—modist{D, B). 


lim inf — log 

T ->00 2t 


IB 


lim sup — log 

r->oo 2 t 


Note that, if we know oq on B, then one can compute the indicator function 


r 



Ve)dx 


from the experimental data and the values of ao on B. 

We just need the following qualitative knowledge: 

(i) the governing equation of the observed wave in the space which has a penetrable obstacle 
takes the form (1.1) and its refractive index is given by (1.2); 

(ii) the governing equation of the observed wave in the space which has no obstacle yet takes 
the form (5.1) with ag in (1.2). 

Summing up, we can say that: one can know the existence of something added to a reference 
space by comparing the “snap shot” u on B with a “reference snap shot” Ve on B even in the 
case when: the reference space has a complicated rough refractive index ao; ao is unknown 
outside B. 


5.2 One-space dimensiona case 

Finally let us describe one non trivial application of the method presented in the proof of 
Theorem 1.1 in the case when the space dimension is one and the background medium is not 
homogeneous. 

We assume that ao is piecewise constant and takes the form 

{ 1, ifx<aor6<x 

(5.2) 

ko, if a < X < b, 
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where —oo<a<b<oo and ko is a positive constant. 

We choose / = xb, where xb denotes the charactersitic function of open interval B = 
]p — e, p + e[ with a fixed p satisfying p + e < a. / is a simple model of the disturbance given at 
t = 0 from the left side of the wall [a, b]. 

We assume that D =]c, d[ with, for simplicity b < c < d < oo. This means that obstace D 
is located behaind the wall ]a, b[ from the observer. 

Define 

p = a - {p + e) + \/ko{b - a) + (c - b). 

The quantity 2(p coincides with the time of flight of the signal which propgates as 


xq = p + e —)• xi = a — X 2 = b —^ x^ = c — X 4 = b —)• X 5 = a —)• xq = p + e, 

where the propagation speed of the signal in ]p + e, a[ and ]b, c[ is 1 , in ]a, b[ is 1 /y/^. 

Note that kowing c is equivalent to knowing p provided b — a which is a thickneess of the 
wall and Xj which is the propagation speed of the wave inside the wall are known. 

Let a be the same as (1.2) with n = 1, a /i G L°°{D) and oq given by (5.2). 

Theorem 5.2. Let T > 2p. We have: 

(i) if h satisfies (A.I), then 


lim e'^'^ 



(a) if h satisfies (A.II), then 


v)dx = — 00 ; 


lim 



v)dx = 00 . 


Moreover, in ease of both (i) and (ii) we have 


T 


lim — log 
^00 2 r ^ 




(5.3) 


The proof is based on the following asymptotic formula of the solution of (1.4) which corre- 
ponds to Lemma 3.1 and the argument of the proof done in that of Theorem 1.1 to the present 
case. 

Lemma 5.1. We have, as r —)■ 00 


2re2^^ / v^dx = l + 0{T-^). 
JD 


(5.4) 


Proof. A direct computation shows that v has the expression 

^(T)e™, if X < p — e 


v{x,t) = { 


5(r)e^^ + C(T)e-^^ + 4, ifp-e<x<p + e 
D(r)e™+ G(r)e-™, ifp + e<x<a 
iL(r)e^^^ + iL(r)e-^™, if a < x < 6 
L(r)e“™, if 6 < X, 
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where 


A{t) = D{t) + 


H{r) 

Kir) 


g-T(p-e) _ e-'r(p+e) 


p-T{p-e) 

B{t) = A{t) -, C{t) = 


2 r 2 ’ 2 r 2 

pT{p+e) _ pT{p-e) 

G{t) = ^ 


2 r 2 


2\/fco ‘^vk{} 

2 \/fco 2\/feo 


(/co - + e-2^^(^-“))(l - 

_ gTbg-TV^(fe-a)g-T(a-p-e) 


pT(p-<^) 

~2^' 


1 1 /\/^-ly (1+ e-2^^(^-“))(l - e-2^") 

\ + 1 / 


Thus, as T —)■ oo 

L(r) = + 0(r-2)). (5.5) 

Since v{x) = L(r)e“'^® for x > b, we have 



_ g- 2 r(d-c))_ 


Then (5.5) yields (5.4). 

□ 

Theorem 5.2 suggests that, if u G //^(R”) is the solution of (1.4) and one knows the leading 
term of ||r’||L 2 (-£)) as r —)■ oo like (5.4) in Lemma 5.1, then one can obtain a formula in three- 
dimensions like (5.3) in Theorem 5.2 instead of estimates (1.5) and (1.6). For the determination 
of the leading term, usually, one has to solve the eikonal equation = Q;o(a;). In some 

restricted cases it is possible to solve the equation, see [18]. The case treated in Theorem 5.2 
is just simplest one of such cases. However, we do not go into such research direction further 
since the solvabilty of the eikonal equation requires some regularity for ao and we are seeking a 
method that works for finding an obstacle embedded in a rough background medium. We think 
that the method presented in this paper is the first one for the purpose. In a forthcoming paper 
we will consider the original problem which is formulated by the Maxwell system. 
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6 Appendix. Proof of Lemma 3.1. 

We make use of the following elementary fact and an iteration process. 

Lemma A. Let A > 0. Given f G L^(R”) let v G be the weak solution of 

{A-X‘^}v + f = 0 inTC. 

Then we have 

ll'W||L2(R,n) < A ^||/||L2(Kn) (^-1) 

and 

l|V^^llL2(Rn) < (2A) ^ ||/||L2(Rn). (^-2) 

Proof of (3.12). 

Rewrite (1.4) as 


{A - + |ao/ + - ao)i'} = 0 inR^. 

Let vi G be the weak solution of 

{A — (Mor)^}ui + aof = 0 inR"^. 


Since vi has the expression 

vi{x)= ao{y)f{y)KMoT{x-y)dy, 

JR3 

we have ui(a;) > 0 . 

Let j = 1, • • •. Given Vj G R^(R’^) let Uj+i G H^(RT) be the weak solution of 
{A - {MoTf}vj+i + {aof + r^(M^ - ao)vj} = 0 inR”. 
Then, Vj+i — Vj for j >2 satisfies 

{(A - {MoTf)}{vj+i - Vj) = -r^(Mo - ao){vj - Vj-i) inR’". 


and from (A.l) we obtain 




and hence, for j = 1 , • 




1-1 


Similary we have 


1112 - 11i||l 2(R«) - l|lll||L2(Rn) 


(A.3) 
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and applying (A.l) to vi on this right-hand side, we obtain 


\v2 - 'yi||L2(R") ^ 1 j ;^II/IIl2(r«). 


Thus we obtain, for j = 1, • • • 

- ^illL 2 (Rr.) < ;^||/||l 2 (R")- 

Moreover, from (A.2) we have 




and similary 


11 ^, ,,, ^ rlMn — mn) .. .. 

||V(uj+i -'Vj)lli,2(R") < --Fi - ^^i-i|lL2(R,«) 


||V(uj+i - ■i;j)||i 2 (R„) < ||/||l 2 (R" 


Therefore, the sequence {uj} in converges to 


V = — Vj) + vi ini/^(R" 

i=i 


(A.4) 


and we have; 

\\v - viWlHh-) < (^1 - II/IIl 2 (r»); 

l|Vu - Vui||i2(R.) < ^ (^1 - ||/||l2(r„). 

Thus taking the limi of (A.3) as j —?> oo, we see that v is the weak solution of (1.4). 

Since V 2 — vi has the expression 

V 2 ix) - vi{x) = T^ (M^ - ao{y))vi{y)KMoT{x - y)dy, 

JR 3 

we have V 2 {x) — vi{x) > 0. For j > 2 we have also 

Vj+i{x) - Vj{x) = tM (Mq - ao{y)){vj{y) - Vj-i{y))KMoT{x - y)dy 

4r3 

and thus by induction we obtain, for all j >2 Uj+i(x) — Vj{x) > 0. Therefore from (A.4) and 
the almost convergence property of a subsequence of {Yl^=i{xj+i — we concude that 

v{x) > vi{x) a.e. X G R"". This completes the proof of (3.12). 

□ 

Proof of (3.13). Let vq G R^(R"') solve 

{A — (mor)^}uo + ao/ = 0 inR”. 


Vq has the expression 

vq{x) = I ao{y)f{y)KmoTix - y)dy 
JR" 
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and thus vq > 0. Then, v has the expression 

V = vo- eo, 

where eo G solves 

Aeo — aoT^€o + t^(o:o — rn^)vQ = 0 inR"'. 

Since T^(ao — m.Q)uo > 0, applying (3.12) to the equation above, we obtain eo > 0. Therefore 

we obtain v < vq. 

□ 

References 

[1] Catapano, I. and Crocco, L., A qualitative inverse scattering method for through-the-wall 
imaging, IEEE Geoscience and remote sensing letters, 7(2010), No.4, October, 685-689. 

[2] Courant, R. and Hilbert, D., Methoden der Mathematischen Physik, Vol. 2 (Berlin: 
Springer), 1937. 

[3] Dautray, R. and Lions, J-L., Mathematical analysis and numerical methods for sciences 
and technology, Evolution problems I, Vol. 5, Springer-Verlag, Berlin, 1992. 

[4] Gilbarg, D. and Trudinger, N. S., Elliptic partial differential equations of second order, 
second.ed. (BerlimSpringer), 1983. 

[5] Hansen, S., The leading singularity of the scattering kernel for a transparent obstacle. Math. 
Ann., 279(1987), 125-139. 

[6] Ikehata, M., Enclosing a polygonal cavity in a two-dimensional bounded domain from 
Cauchy data. Inverse Problems, 15(1999), 1231-1241. 

[7] Ikehata, M., Reconstruction of the support function for inclusion from boundary measure¬ 
ments, J. Inv. Ill-Posed Problems, 8(2000), 367-378. 

[8] Ikehata, M., Extracting discontinuity in a heat conductive body. One-space dimensional 
case. Applicable Analysis, 86(2007), no. 8, 963-1005. 

[9] Ikehata, M., The enclosure method for inverse obstacle scattering problems with dynamical 
data over a finite time interval. Inverse Problems, 26(2010) 055010(20pp). 

[10] Ikehata, M., The enclosure method for inverse obstacle scattering problems with dynamical 
data over a finite time interval; H. Obstacles with a dissipative boundary or finite refractive 
index and back-scattering data. Inverse Problems, 28(2012) 045010(29pp). 

[11] Ikehata, M., The enclosure method for inverse obstacle scattering problems with dynamical 
data over a finite time interval: HI. Sound-soft obstacle and bistatic data. Inverse Problems, 
29(2013) 085013(35pp). 

[12] Ikehata, M., Extracting the geometry of an obstacle and a zeroth-order coefficient of a 
boundary condition via the enclosure method using a single reflected wave over a finite 
time interval. Inverse Problems, 30(2014) 045011(24pp). 


19 



[13] Ikehata, M., New development of the enclosure method for inverse obstacle scattering, to 
appear as a chapter in Inverse Problems and Computational Mechanies (eds. L. Marin, L. 
Munteanu, V. Chiroiu), Vol. 2, Editura Academiei, Bucharest, Romania. 

[14] Ikehata, M. and Itou, H., On reconstruction of a cavity in a linearized viscoelastic body 
from infinitely many transient boundary data. Inverse Problems, 28(2012) 125003(19pp). 

[15] Majda, A. and Taylor, M., Inverse scattering problems for transparent obstacles. Comm, 
in partial differential equations, 2(4)(1977), 395-438. 

[16] Petkov, V., Propagation of sigularities and inverse scattering problem for transparent ob¬ 
stacles, J. Math. Pures Appl., 60(1982), 65-90. 

[17] Petkov, V., Inverse scattering problem for transparent obstacles. Math. Proc. Camb. Phil. 
Soc., 92(1982), 361-367. 

[18] Uhlmann, G., The Cauchy data and scattering relation, in Geometric methods in inverse 
problems and PDE control IMA vol. 137(New York: Springer), pp.263-288, 2004. 

[19] Zetik, R., Crabbe, S., Krajnak, J., Peyerl, P., Sachs, J., et ah. Detection and lo¬ 
calization of persons behind obstacles using M-sequence through-the-wall radar, Proc. 
SPIE 6201, Sensors, and Command, Control, Communications, and Intelligence (C3I) 
Technologies for Homeland Security and Homeland Defense V, 620101 (May 10, 2006); 
doi:10.1117/12.667989. 


e-mail address 

ikehat a@ amat h. hiroshima-u. ac. j p 


20 



